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Abstract: The Peccei-Quinn (PQ) mechanism provides an elegant solution to the strong 
CP problem. However astrophysical constraints on axions require the PQ breaking scale to 
be far higher than the electroweak scale. In supersymmetric models the PQ symmetry can 
be broken at an acceptable scale if the effective potential for the pseudo-modulus in the 
axion multiplet develops a minimum at large enough field values. In this work we classify 
systematically hadronic axion models in the context of gauge mediation and study their 
effective potentials at one loop. We find that some models generate a PQ scale comparable 
to the messenger scale. Our result may prove useful for constructing full realistic models 
of gauge mediation that address the strong CP problem. We also comment briefly on the 
cosmological aspects related to saxion and axino, and on the quality of the PQ symmetry. 
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1 Introduction 

The Peccei-Quinn (PQ) mechanism [1, 2] is one of the most plausible solutions to the strong 
CP problem [3]. It involves an anomalous U(1)pq symmetry that is spontaneously broken. 
The corresponding pseudo-goldstone boson is the axion [4, 5] . Astrophysical constraints [6] 
on axions imply that the PQ breaking scale has to be very high {fa > 10® GeV). This 
introduces a large hierarchy between the PQ and the electroweak scale which calls for a 
theoretical explanation. One way to address the issue is in the context of Supersymmetry 
(SUSY), using the SUSY breaking effects to generate the PQ scale dynamically. 

This possibility was first considered by two different groups [7, 8]. They proposed a 
supersymmetric hadronic axion model composed of three sectors: a hidden sector where 
SUSY is broken, a PQ sector defined by the superpotential 

VPpQ = , (1.1) 

and a messenger sector 

Wm = ■ ( 1 - 2 ) 

The chiral superfields S, 4>p and are charged under the U(1)pq, with S the gauge-singlet 
axion multiplet, and <hp, 4>p in the 3 and 3 of the color gauge group SU(3)c- The superfield 
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X has a SUSY breaking vacuum expectation value (VEV), {X) = M + 6‘^F, and 
are the messenger superhelds of gauge mediation, also in the 3 and 3 of SU(3)c'- At two 
loops gauge mediation generates SUSY breaking masses for the PQ squarks and In 
turn SUSY breaking effects are transmitted to the scalar S field via one-loop diagrams of 
and "hp. This generates an effective potential V'eff('S'), effectively at three loops, such 
that < 0. The effects of gravity mediation, parametrized as 

UO 

Ug,~m2/2|5|2, (1.3) 

with m ^/2 the gravitino mass, then stabilize the potential for S at large field values. This 
results in a large PQ breaking scale, fa = {S). 

A similar setup was reexamined more recently in the context of gauge mediation. 
Some of the present authors considered a model in which the PQ messengers <hp, $p and 
the regular messengers <I>m, mixed [9] and concluded that at one loop it did not stabilize 
the PQ scale away from the origin of the field S. The same work [9] considered other ways 
of generating the PQ scale dynamically, but they all involved models with additional gauge 
interactions that were somewhat complicated. In Ref. [10] the authors considered a slightly 
different model with messenger mixing and showed that it led to the stabilization of the 
PQ scale at two loops. ^ 

In this work we revisit the dynamical generation of the PQ scale in the context of 
gauge mediation and we ask if it is possible to achieve it in a simpler way. To this end we 
consider models with an i?-symmetry and a PQ symmetry, containing the fields S, X and 
an arbitrary number of messengers. We classify these models systematically according to 
the charge assignments of the fields, and point out possible virtues and disadvantages for 
each class. We study in detail concrete examples in which the one-loop effective potential 
generates a VEV, {S) ~ M, with M the messenger scale. Taking a high enough messenger 
scale, M > 10® GeV, we then have an acceptable PQ breaking scale, fa ~ M. This is 
significantly simpler than the scenarios we described above. First, the analysis at one loop 
is sufficient. Second, the PQ scale is stabilized solely within the context of gauge mediation, 
without the need of intervention of gravity mediation. 

The paper is organized as follows. In section 2 we introduce the models and cate¬ 
gorize them into four distinct classes. In Section 3 we study a couple of examples and 
comment on the vacuum stability and on cosmological constraints. In Section 4 we discuss 
how one could address the issue of the axion quality in these models by invoking discrete 
symmetries. In section 5 we summarize our main results and discuss possible implications 
for phenomenological model building. We include two appendices. In the first we report 
many details on models with two sets of messengers. In the second we study the anomalies 
related to the discrete gauge symmetries, invoked to address the axion quality problem. 

2 Classification of models 

We consider models defined by the superpotential 

W = M^j(X, + Wn(X) , (2.1) 

^See also Refs. [11, 12] for different mechanisms of stabilization of the PQ scale in gauge mediation. 
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Table 1: R and PQ charges of the fields in the model. 


where 

Mij{X,S) = XXij+mij + S5ij , i,j = 1,...,N. (2.2) 

There are N messenger fields, transforming under 5 and 5 of SU(5) respectively, and 

two gauge-singlet fields, X and S. The matrices A, m and 5 are N x N. We assume a global 
/^-symmetry, U{1)r, and a PQ symmetry, U(1)pq, under which the fields are charged as 
shown in Table 1. The choice of keeping X neutral under U(1)pq is convenient as it keeps 
the SUSY breaking and PQ breaking sectors separate. Note also that our models, as most 
models of gauge mediation, have an extra U{l)v global symmetry, under which the d>’s 
and the <h’s transform with opposite phases. The quantum number associated with this 
symmetry is the messengers’ number.^ 

The term W/{(Y) in Eq. (2.1) stands for a hidden sector superpotential, necessary to 
generate a SUSY breaking VEV for X. A minimal choice is Wr(Y) = FX but for the 
time being we leave Wr{X) unspecified. Consistent with our choice that X is uncharged 
under the PQ symmetry, we assume that all the fields that appear in Wr{X) are singlets 
under U(1)pq. 

Our notation, as well as the reasoning we outline in the rest of this section, are inspired 
by the models of (extra) ordinary gauge mediation [14] where the superpotential has the 
form of Eq. (2.1) with S absent. The superpotential in Eq. (2.1) must have R charge 2 and 
PQ charge zero, resulting in the following selection rules 


\ij / 0 

5ij / 0 


only if fi + rj = 0 
only if fj + rj = 2 
only if Vi + rj = 2 — rs 


and Pi+Pj = 0 , 
and Pi+Pj = 0 , 
and Pi + pj = -ps ■ 


As a consequence the determinant of the matrix A4 is a monomial in X and S: 


(2.3) 

(2.4) 

(2.5) 


detM = Y",S'?G(A,m,5), 


( 2 . 6 ) 


where G(A, m, 6) is some function of the couplings and 

N 

n = 

i=l 




PS 


q = 


1 ^ 

-^{Pi+Pi)- 

ley < ^ 


P^t=i 


(2.7) 

( 2 . 8 ) 


the messengers’ number were conserved the lightest component of the messengers would be stable and 
could overdose the universe. We assume that U (l)v is broken in another sector of the theory. A coupling of 
one messenger to two matter fields in the superpotential of the visible sector, for instance, would suffice [13]. 
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Here n and q are integers satisfying 0 < n < iV and 0 < q < N. The proof of the identity 
in Eq. (2.6) can be done in analogy with that given in Ref. [14]. 

Our aim is to look for models that generate a large VEV for the held S, so that the PQ 
symmetry is spontaneously broken. Note that this implies that the PQ charge of S must 
be strictly nonzero, ps ^ 0. As 5 is a hat direction of the classical potential we have to 
compute the one-loop effective potential and check if it stabilizes S away from the origin. 
The identity of Eq. (2.6) leads to a classihcation scheme with four qualitatively distinct 
types of models. 

• Type A: det m / 0, det A = det <5 = 0. 

Here we can go to a basis where m is diagonal, which implies Vi+Vi = 2 and Pi+Pi = 0. 
It follows that 

n = q = 0 . (2.9) 

The messengers are stable around X = 0 and 5 = 0, but some can become tachyonic 
at large X and S. In these models the gaugino masses vanish to leading order in F. 
We will not study any models of this kind in the remainder of the paper. 

• Type B: det A / 0, det m = det J = 0. 

Here we can go to a basis where A is diagonal, f j -|- r* = 0 and Pi + Pi = 0. Thus we 
have 

n = N and q = 0. (2-10) 

At large X all messengers have masses of order AX. As X approaches the origin, 
detm = 0 implies that some messengers have 0{m) masses, while others are light 
with masses going to zero as some power of X. Eventually the latter become tachyonic 
as X gets closer to zero. There can however be local minima of the potential for finite 
X where all the messengers are massive. We study a simple model of this form in 
Appendix A. 

• Type C: det 5 ^ Q, det A = det m = 0. 

Here we can go to a basis where 5 is diagonal, ri + ri = 2 — rs and Pi+Pi = —ps- It 
follows that 

n = 0 and q = N. (2-11) 

With a reasoning analogous to the one for Type B models, we deduce that in Type C 
models some messengers become tachyonic as S goes to zero. As n = 0, the gauginos 
are massless at leading order in the SUSY breaking parameter F. In the next section 
we study in detail a model belonging to this category for which S is stabilized away 
from the origin. 

• Type D: det A = det m = det <5 = 0. 

In this category we have 


0 < n < N and 0 < q < N. 


( 2 . 12 ) 
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The messenger sector combines features of the previous categories and there are no 
tachyons in a region Xmin < |-^| < -^max, *S'min < |‘S'| < S'max- We will see that models 
of this kind can develop a minimum with {X) ^ 0 and {S) 7 ^ 0 at one loop. 


3 Concrete examples 

3.1 A model with two sets of messengers 

In this section we consider a model of Type C with two sets of messengers {N = 2) that 
generates a VEV for the field S of order the messenger scale M. Taking M > 10® GeV 
we then have an acceptable PQ breaking scale. Many details of the analysis are found 
in Appendix A where we look in general at models with N = 2 and find other options 
(of Type B and D) with less appealing features. For the sake of definiteness we specify 
which we take to be of the simple form proposed in Ref. [15]. The superpotential 
of the model is then given by: 


W = W^q + Wr-, (3.1) 

+ SXs , (3.2) 

Wr = X{X(fi(p-i + F) + mnp-npz + ■ (3-3) 

Here the (/? fields are gauge- and PQ-singlets, with the subscript denoting their R charge. 
The PQ and R symmetries forbid additional renormalizable terms in Eq. (3.1). We have 
taken the two couplings of S to the messengers to be equal. With this choice the model has 
a messenger parity symmetry [16] that has the virtue of forbidding dangerous D-terms [17], 
which would otherwise lead to tachyonic sfermions. 

At tree level there is a pseudo moduli space of vacua on which Lpi = 0, = 0, $* = 0 

with X and S arbitrary. As described in the previous section some of the messengers are 
tachionic for small S. Indeed for small S the potential rolls down to a moduli space of 
supersymmetric vacua with = 0, 5 = 0, X = 0 and on which the gauge invariant 
combinations are subject to the constraints: 

F ~ ~ 

$ 1^1 = -—, $i^> 2 -f ^> 2 $i = 0 . (3.4) 

Moreover on the pseudo moduli space some of the (pi become tachionic at large X. In 
this case the potential rolls down along a runaway direction on which = 0. This 

runaway is parametrized by —)• 00 and [15] 


V >?F ) 




( Frriiips X 3 

V Am2 J 


/ F‘^m2 \ 3 

yX'^miipsJ 


(3.5) 


We are interested in establishing if the one loop effective potential on the pseudo moduli 
space has a local minimum in X and S such that no held is tachionic. The Coleman- 
Weinberg formula for the potential at one loop is: 


p(i) 


Gdvr^ 


STr 






(3.6) 
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where STr stands for supertrace, is shorthand for the squared mass matrices of the 
scalar and fermion components of the superfields and A is the cutoff scale. We can compute 
the one-loop effective potential S'), as a function of both X and S', at all orders 

in the SUSY breaking parameter F. Given the form of W in Eq. (3.1) the squared mass 
matrices of scalars and fermions are block diagonal and the resulting one-loop potential 
can be written as the sum of two terms 


U(i)(Y,5) = Uj(^)(Y,5) + V^^\X ), (3.7) 

where V^\x) was computed in Ref. [15]. For some range of the parameters one obtains a 
minimum at {X) = M for which all the (fi's are non tachyonic. The IpQ {X, S) contribution 
to the one loop potential will not destabilize this minimum provided that 


dvi]l(X,S) 

dX 


«(Y) 


d^v^^\x) 

dX^ 


(3.8) 


This condition can be satisfied by taking \x sufficiently small, < 0.1 A. 

In Appendix A we present a detailed analysis of Upq(Y, S). In summary there is a 
local, PQ breaking minimum at 




(3.9) 


For S < Stac, with Atac = some messengers become tachyonic and the system rolls 

down classically to the SUSY vacuum. For Amin to lie in the tachyon-free region we need 
to satisfy 

F<e-^XxM^. (3.10) 

We want to make sure that our metastable vacuum is long lived, i.e. that the tunnelling 
rate from Amin to A = 0 is low^. A detailed study of this rate is beyond the scope of this 
work. It is sufficient to note that the potential barrier height does not vary strongly with 
As while its width, Amin — Atao is proportional to Aj^. Hence the tunneling rate can be 
made small by lowering As. 

At the metastable minimum of UA)[x. A), with (A) ~ {X) = M, we have the following 
spectrum: 


1. The axion and the R-axion are massless.^ 

2. The saxion has a mass 

/ A2 F 

In terms of loop counting this is larger than that of the MSSM particles in gauge 
mediation. However vacuum stability requires small As, and as a result the saxion is 
likely lighter than most of the MSSM particles. 

®The tunneling rate to the runaway direction where ipa —>■ oo can be easily made very small [15]. 

^The axion then acquires a small mass due to the U(1)pq anomaly, while the J?-axion can acquire a 
mass from explicit it-symmetry breaking in supergravity [18]. 
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3. The axino, the fermionic component of the superfield S, has a mass 


rrin 




F 


Ihvr^ M ’ 


so it is lighter than the saxion. 


(3.12) 


4. The ii-saxion has a mass 


rriRs ~ 


A2 F 
IGtt^ M ’ 


(3.13) 


and is typically heavier than the saxion, as the coupling A that appears in Wr does 
not need to be as small as A*. 


5. The gravitino is light, as usual in gauge mediation 


m3/2 ~ 


F 


where Mpi is the reduced Planck mass. 


(3.14) 


The saxion and i2-saxion are the two pseudomoduli in the model and we need to make 
sure that they do not pose cosmological issues [19]. Typically such issues are more severe 
the lighter the pseudomodulus, so we concentrate on the saxion here, as mg < niRg. The 
main decay modes of the saxion are (i) into two axions, (ii) into an axino and a gravitino, 
and (iii) into two gravitinos.® The decays (ii) and (iii) can be understood from the effective 
operator 

^ J d^e{X^X){S^S). (3.15) 

The decay rate in each case is given approximately by 


r 


s 


r\j 


1 rrig 

IGtt m2 


~ 10 


-25 


GeV 


/ rug \3/l0^2QgY\ 

Vl GevJ V M ) 


(3.16) 


The saxion starts oscillating about its minimum when the Hubble rate, H r\j r2/Mpi, is 
comparable to its mass, H ~ mg. At that time it has an energy density of order 
and constitutes a fraction /Mp^ of the total energy density. From then on it behaves 
like matter, so its fraction of the energy density grows with the scale factor. It decays when 
if ~ Ts, that it is at a temperature 


t: 


idee 


3/2 

mg' Mp] 

10 M 


(3.17) 


Requiring that it decays safely before Big Bang Nucleosynthesis (BBN), > 0.1 GeV, 
amounts to a lower bound on XgF: 


XgF > 10^5 GeV^ 


M 

1012 GeV 



(3.18) 


®The saxion can also decay into pairs of MSSM particles or pairs of gauge bosons, but these are further 
suppressed and subdominant compared to (i), (ii), and (iii) in the text. 
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where we used Eq. (3.11). As long as M < 10^^ GeV, the bound of Eq. (3.18) also 
guarantees that the saxion decays before it comes to dominate the energy density. This 
is helpful in relieving possible constraints from extra radiation at BBN [20], as the decay 
products are relativistic. 

The axino in this model is much heavier than the gravitino and its main decay mode 
is into a gravitino and an axion, which are both dark matter candidates. Avoiding over¬ 
closure of the universe with too much gravitino dark matter then results in a bound on the 
reheating temperature, which should be lower than about 10^ GeV [ 21 ]. 

3.2 Models with more messengers 

The model considered in Eq. (3.1) is of Type G and has the deficiency that gaugino masses 
are suppressed, as they are not generated at one loop. One way to fix this is to introduce 
more messengers. For example we can replace Wr in Eq. (3.1) with 

Wr = AA($3<1>3 + I> 4 <h 4 ) + XF + ml> 3 <l >4 . (3.19) 

Under the i?-symmetry the new messengers have charges r 4 = r 3 = 1 and f 3 = r 4 = — 1. 
The model specified by this Wr plus WpQ [which we take of the same form as in Eq. (3.2)] 
is of Type D. The superpotential (3.19), which was considered in Ref. [14], breaks both 
SUSY and the R-symmetry spontaneously for a large range of parameters. This continues 
to be true when added to the PQ sector as long as the couplings in VU^q are not too large.® 
This construction is rather ad hoc but serves as an example of a messenger sector that 
breaks SUSY, the R-symmetry and the PQ symmetry at one loop. 


4 The axion quality 

The PQ mechanism provides a solution to the strong GP problem as long as the U(1)pq is 
violated almost exclusively by the QGD anomaly. One must ensure that higher-dimension 
operators which violate the symmetry explicitly [22, 23] are suppressed. In our framework 
such operators can appear in the superpotential, taking the form 

S 

) where a G N®, 6 G N"*", (4.1) 

and Mp is the Planck mass. These terms result in a contribution to the scalar potential 
that has to be small: 



Here we have taken {X) ~ (S) = fa — M. The above inequality must be satisfied if the 
axion is to provide the solution to the strong GP problem [24]. If fa — 10^® GeV and 
\F\ ~ 10^® GeV^, then this corresponds to the requirement that operators with a -|- 6 < 8 
must be forbidden in the superpotential. 

®A11 possible couplings between the sets of messengers appearing in WpQ and the messengers in Wr 
need also be small. 


SW D XMi 


{-] 

\Mp 
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Additionally one should worry about non-renomalizable operators in the Kahler po¬ 
tential. These are of the form 


5K D 


d^ed^exx^ aGN°,6GN+, (4.3) 


which leads to the following contributions to the scalar potential 

a-\-b 


FF^ 

\Mp ) 


< 10"^^ GeV^. 


(4.4) 


Taking once again fa — 10^*^ GeV and 1^1 ~ lO^^GeV^, we must forbid operators with 
a + b < 5. This is a less stringent constraint on a and b compared to contributions from 
the non-renormalizable superpotential operators of Eq. (4.1). 

One possibility to address this problem is to impose discrete gauge symmetries [25-27]. 
We consider imposing a discrete PQ symmetry, and a discrete /^-symmetry, Z^, such 
that the global U(l)ij and U(1)pq described in the previous sections arise as accidental 
symmetries of the Lagrangian. The dangerous operators of Eq. (4.1) and Eq. (4.3) then have 
to respect Z^^ x Z^ and we have to find how large N and M need to be in order to satisfy 
Eq. (4.2). The constraints arising from the requirement that these discrete symmetries are 
anomaly free are discussed in Appendix B. The analysis is model dependent as it involves 
details of the visible sector (or other sectors) of the theory. As an example, for the model 
defined by Eq. (3.1), some of the possibilities are: 


1. We can neglect anomaly constraints in our hidden sector model [Eq. (3.1)] under the 
assumption that extra matter charged under SU(5), Z^^ and Z^ is introduced to 
cancel anomalies [28]. The minimal suitable discrete symmetry is then Zg^ x Z^, 
with the discrete R charge of S chosen as = 1. 

2. If there is no additional matter charged under Z^^ and SU(5) we can still neglect 
anomaly constraints on the Z^ symmetry by assuming that the R charges in the 
visible sector are chosen to cancel any resulting anomaly. The minimal suitable 
discrete symmetry is then Z^^ x Z^, with =2. 

3. We can assume that the visible sector is anomaly free by itself. The minimal com¬ 
patible discrete symmetry is then Z^*^ x Z{\, with qg = 5. 


5 Summary and discussion 

The feasibility of the axion solution to the strong GP problem requires a PQ breaking scale, 
fa, much above the electroweak scale. It would be desirable to explain this hierarchy with a 
mechanism that generates fa dynamically. The main result of this paper is the construction 
of simple models that tie fa to the messenger scale M of gauge mediation. This is achieved 
with a one-loop analysis where it is sufficient to compute the Coleman-Weinberg potential, 
induced by SUSY breaking effects, for the scalar component of the axion superfield S. We 
find some examples in which such a potential leads to {S) ~ M. This implies that the 
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messenger scale, M ~ fa, has to be larger than 10® GeV, which puts it on the higher end 
of the range typically considered in gauge mediation. 

Our models possess an i?-symmetry and a PQ symmetry, and can be classified in a way 
similar to Ref. [14]. We find four distinct classes but study in detail only a few examples 
(some of which are in Appendix A). There is a lot left to explore. For instance, we have not 
studied any model of Type A. They suffer from the issue of suppressed gaugino masses, but 
perhaps one can address that problem in another sector of the theory. Models of Type D 
do not have such an issue and are possibly the most interesting to further investigate. We 
have mentioned only one example of a successful Type D model with four sets of messengers 
in Section 3. One can be more clever and find additional examples where the same sector 
spontaneously breaks the i?-symmetry (thus breaking SUSY) and the PQ symmetry. This 
would specify the entire “hidden” sector, which could then be connected to the visible 
sector in the gauge mediation framework. The result would be a full calculable model for 
which one could study the phenomenological implications. 

We have also addressed the problem of the axion quality in Section 4. We do not 
provide any new insight into this issue, and show that relatively large discrete symmetries, 
thus somewhat unattractive, are needed to ensure the high quality of the PQ symmetry. 
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A Classification of models with two sets of messengers 

The simplest subcategory of the models we introduced in Section 2 consists of models with 
only two sets of messengers. As the goal is to generate a VEV for S via the SUSY breaking 
effects encoded by {X), we need at least one nonzero entry for Ajj and one nonzero entry 
for 5ij. With N = 2 this immediately implies m = 0, which means that models in this 
subcategory can only be of Type B, C or D. In this Appendix we show that with this field 
content there is one model for each Type (except A) which is representative of any model 
one could write. 

We assume a canonical Kahler potential for all the fields, and we write the superpo- 
tential as 

IT = 1Ur(A) + lUpQ . (A.l) 

We start from the most general Wpq consistent with the SU(5) and R symmetries: 

lUpQ = A + Xx22^2^2 + (A.2) 

-h S (^A^„l>i$i Xs22^2^2 + A^i2l>i$2 + A^2iI>24>i) . (A.3) 


10 


We do not write the term XF here, also consistent with the symmetries, but we assume it 
is included in Wr{X). We take all the couplings A’s to be real. We impose a global U(1)pq 
under which S and the messengers are charged but X is neutral: 


(A.4) 

We then have the following equations for the PQ charges: 

Pi+Pi = 0 Pi+Pi = -Ps (A.5) 

P2+P2 = 0 P2+P2 = -Ps (A.6) 

Pi + P2 = 0 p^+p2 = -ps (A.7) 

P 2 +Pi = 0 P 2 +Pi = -Ps- (A.8) 


Clearly the two equations on each line are mutually exclusive. If we satisfy three equations 
of either column, the remaining one of the same column is also automatically satisfied, 
which would allow only operators with either X or S'. We need at least one operator with 
X and one with S in order to generate the PQ scale through the SUSY breaking effects. 
Thus we can satisfy at most two equations of either column. 

Suppose we keep the hrst two equations on the right column. Then it is easy to check 
that we can only satisfy either the last or the second to last equation on the left column. 
If we keep the second to last we have 

XA,,,$i$ 2 + s (a,,,I>i4>i + A,,2^>2^'2) . (A.9) 

Here we can relabel the unbarred fields d>i -H- d> 2 , so that the above is equivalent to 

W^Q = XA^,,l>i$i + S (a,,,I>i$2 + A,,,I>2$i) . (A.IO) 


Keeping the equations on the right of (A.7) and (A.8) and the one on the left of (A.5) is 
equivalent to keeping (A.5) and (A.6) on the right and (A.7) on the left as we did above. 
Had we chosen instead (A.5) and (A.6) on the right and (A.8) on the left, we could have 
relabelled -H- 4>2 and would have ended up again with the superpotential in Eq. (A.IO). 
This defines N = 2 models of Type C, with X coupled diagonally to only one set of 
messengers and S coupled off-diagonally to both sets of messengers. 

With a procedure analogous to the one just outlined, we hnd N = 2 models of Type 


B: 


W^Q = X + Ax22 ^2$2) + SXs,2^1^2 ■ 


(A.ll) 


Next, we can keep the left of (A.5) and (A.7) and the right of (A.6) and (A.8), which 
leads to 

X + A,,i2$i$2 ) + S (^A,22«l2^2 + Xs2i^2<^i) • (A.12) 

Here one of the operators is redundant. We have pi = P 2 , so we can use the field redefinition 
to get rid of the second operator involving X. After relabelling fields and couplings we 
have 

W^Q = XA,,,,I>1$1 + 5 (a,22 ^2$2 + A.2 i^2$i) . (A.13) 
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We could have chosen a different field basis. For example we could have used the field 
redefinition to get rid of one of the operators with S instead. Of course the resulting 
physics would be equivalent. 

Keeping the left of (A.5) and (A. 8 ) and the right of (A. 6 ) and (A.7) we have 


A + 5 (a,,2 -^2^2 + A,,2 «Ii$2 ) . (A.14) 

Here pi = p 2 and we can redefine the fields <hi and $2 to get rid of one of the operators. 
Again we get rid of the second one and we write 

+ 5 (a,22 ^>2^2 + • (A.15) 

This is equivalent to eq. (A. 13) under the exchange <l>j -H- 4>j. 

To summarize, we have seen that different choices of the PQ charges lead to: 

= A (a,,,I>i4>i + X.,,^ 2 ^ 2 ) + SXs,,^i^2 (A.16) 

= AA,,,,l>i$i + 5 (a,,,4>i$2 + A.,,4>2$i) (A.17) 

= AA,,,,l>i$i + 5 (a,,,I>2$2 + A,,,I>2 $i) . (A.18) 

The scalar potential obtained from Hpq has a pseudo-moduli space given by 

<l>i = = 0 , A and S arbitrary. (A. 19) 


The flat direction along A is lifted by the one-loop dynamics of 1 Tr(A). This in turn leads 
to lifting the fiat direction along S. We can check if for any of our three models in this 
Appendix the one-loop potential in the S'-field direction has a minimum which is not at 
the origin. That would lead to the spontaneous breaking of U(1)pq. Also we would like 
the 5-field VEV to be large in order to generate an acceptable PQ scale. 


A.I One-loop potentials at order 

Under the assumption that F <C M^, instead of using Eq. (3.6) we can study the one-loop 
potential up to order O(F^) with a simpler formula [29] 


V-m = -K.J = ^ y TV (m'M log ^) . 

Here the mass matrix A4 is defined by 

kUpQ = (^1 ^2)-X4 



(A.20) 


(A.21) 


M = A4(A, 5) is a function of the fields A and 5. The computation of the potential of 
Eq. (A.20) proceeds as follows: 

1. compute the eigenvalues of and take the trace in Eq. (A.20); 
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2. substitute X with its VEV {X) = M + 6^F and expand the integrand for small F 
{F < m2); 

3. perform the integration. 

This method of calculating the one-loop potential allows us to examine each one of our 
models without specifying the explicit form of Wji- 

For the rest of this section X and S denote the scalar components of the corresponding 
superfields. Once we determine VpJ{S) we want to check whether it produces a large {S). 
The expressions we obtain from eq. (A.20) with arbitrary couplings A’s are complicated 
enough not to allow us to calculate analytically the extrema of the potential. Therefore we 
proceed in two steps: 

1. we expand 1 /^ 2 ^(S') for l^l <C M and we check if we have a local maximum or minimum 
at S' = 0; 

2. we expand 1 /^ 2 ^ (S) for |S| 3> M. In all the models we find that the potential increases 
at large values of S. 

If there is a maximum at S = 0 we conclude there must be a minimum at S 7 ^ 0. If there 
is a minimum at S = 0 it could be local, and in principle one could have a global minimum 
at S / 0, although this never turns out to be the case for the models we study. 


A.1.1 Type B Model 

Eq. (A. 11) yields a mass matrix of the form 


\ 0 Xx22^ , 


(A.22) 


This model has been considered in Ref. [9]. The one-loop potential at order F'^ i 


IS 


y(i) _ 

Vp2 - 


327r2 


2C 

IFB 


+ ( + Xl^ - In [a - 


+ ( A^,, + Xl^^ + ^6^3/2) In 


m 2 


A + R^/ 2 ^ 


where 


A — Aj,^^ -I- Xx22 + 


ALI^I' 


M 2 


B = A^-4XlX^^, 


C = M^Xl, - XlJ^Xl^ + a2^J + A2^JS|2(A4^^ + 6A2^^a2^, + Xtj, 


(A.23) 


(A.24) 

(A.25) 

(A.26) 


D = {Xl^, - Xl22fM^ + XIJSWXI^^ + Xlj + A^^^|S|"m2(3A^^^ - 2A2^^a2^^ + 3Xtj 


+ 3A^^2|S| M {Xx^^+Xx^2 


)(AL - A^22)^- 


(A.27) 
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Type B Model 



Figure 1: A plot of the full Kahler effective potential (blue, Eq. (A.23)) and the expansion 
of this potential around the origin (brown, Eq. (A.28)) for the Type B model. We set 
■^xii = = 1, Ax 22 = 0-9, A = M, and F = 


For small IS"! 


hhisi« M) = 




^2 x2 . 2,1 A" /A2 M" 

A„,, + A™,., + In ( ——^— I + In 


+ 


327r2 


A 2 




AS -A 2 

V'a;ii ''X22/ 


A 2 


\s\\^f\s\^ 


M4 


M 6 


(A.28) 


the first three derivatives are zero at IS"! =0 and 

dS^dS*'^ 87r2A2^^M4^^^^ 

Here we have defined y = A 3 ; 22 /Aa;i 2 and 

— 2y‘^ Iny^ — 1 


(A.29) 


(A.30) 


The function c{y) is positive for any y, thus we have a local minimum at the origin. 

At large IS"! the potential grows logarithmically 

V'“ (|S| » M) = ^ In ° (j^) ■ 

We can check graphically in Fig. 1 that the minimum at the origin is the global one. Thus 
this model does not seem to generate a VEV for the field S, as already stated in Ref. [9]. 
Such a conclusion, however, is based on the analysis at order F^. 
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Type B Model - Expansion in F Type B Model - Expansion in F 




Figure 2: The effective potential for model II. The green curve corresponds to the full 
effective potential, i.e. all orders in F, the blue is up to order the red is up to order 
F^. On the left-hand plot the green curve lies almost directly behind the red curve. On 
the right-hand side we plot the same potential over a smaller domain: the minimum and 
the difference between the full potential and the F^ expansion become apparent. Here a 
relatively large value of F is chosen, namely F = 0.2M^, and we set A = M. 


Let us see what happens when we consider higher orders in F. To perform the analysis, 
as we did in Section 3, we first have to specify Wr{X). Again, we use the one proposed 
in Ref. [15]. Then we can use Eq. (3.6) to compute the effective potential. For the sake of 
simplicity we take all the couplings to be equal = Xx 22 = "^*12 = "^^*21 = 1- 

Shown in Fig. 2 is the full one-loop effective potential and its expansion to 0{F^). 
We see that there is indeed a minimum not located at the origin. This occurs only when 
including terms beyond F ^, and was missed in the analysis of Ref. [9] . To find the location 
of the minimum, one can series expand the full one-loop potential both in S and in F 
yielding the following expression 


i?2 


/M2 

SVT® 



to 

2|5|^ 

12M4 1 10M4 

M2 


+ 0 


1 ^ 

M6 


(A.32) 


Setting the first derivative with respect to S to zero we find 


(|.?l) ^ ^ (A.33) 

To obtain an acceptably high PQ scale in this model we need F/M > 10® GeV, which also 
leads to a very heavy supersymmetric particle spectrum. 
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A.1.2 Type C Model 

This is the model we studied in detail in Section 3. From Eq. (A. 10) the mass matrix is 


C _ ( As^2*S' 


= 


\s2lS 0 


(A.34) 


The one-loop potential at order F'^ 


IS 


y 2 p2 

_ 3=11 

327r2 


di) 




B 


^3/2 


In 


V 2 A 2 


A- 


^ 1 / 2 ^ 


+ 


+ A)(B- A^^^A + 


B3/2 


V 2 A 2 


A+ 51 / 2 ^ 


where 


A = 


I^P(AL + AL) , .2 


A/2 
4|5|"A2 A 2 


+ A 


Xii’ 


JJ -I - I - 'S12 'S21 I 

®-Mi-+ 


-I- Ag^J 

Af 2 


+ a: 


Xll 


(A.35) 

(A.36) 

(A.37) 


At small S we find 
di) 


/?2 \2 

r'S'dSi « M) = 


1 + In ( Aj.j^ -TTT 


A 2 ) 


+ 




F 2 |5|4 
327r2A2^^ W 

aLa^ 2 i |5|" 


Asi2 + + As2i 


•§12 -§21 _ 

Al, M4 


+ 0 


M 6 


. (A.38) 


Note that there is no term proportional to \S\'^. The stability of the origin is dictated by 
the sign of the fourth derivative 


94y^2^(|5| <CM) 
95255*2 


i ?2 




Xi,,+20Xl^Xl^+Xi^+2Xl^Xljn 


aLa« 2 i 


■§12 -§21 _ 

Al, M4 


(A.39) 


This is negative for |5|/Af —)• 0, as the logarithmic term dominates over the other positive 
contributions. Thus we have a local maximum at the origin. 

At large |5| we have 


E^' 2 ^ (|5| » M) 


2\2 


F^X 


167r2 


In 




(A.40) 


so the potential increases monotonically. Therefore this model has a minimum away from 
the origin. From Eq. (A.38), setting = 0 , we find 


(|5|) ~ ±M 


'xii 


(A. 


A 


exp 




Afi2 + At2i \ 

2 ALA^ 2 i )\ 


(A.41) 
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Type C Model 



S/M 

Figure 3: A plot of the full Kahler effective potential (blue, Eq. (A.35)) and its expansion 
around the origin (brown, Eq. (A.38)) for the model of type C. We set = -^S 2 i = 

1, A = M and F = 


The equality is not exact as we used the expansion for jS*! <C M, but this result is useful 
for some analytic guidance. Setting A^^j = -^*21 = A^u = A^, and defining 


X = 


(A.42) 


we have 


(|5|) ~ ±x 


(A.43) 


We plot in Fig. 3 the full one-loop potential in blue and its expansion around the 
origin, Eq. (A.38), in brown. We see that the actual minimum is at a slightly lower value 
than that of Eq. (A.41). Its location can be adjusted by playing with the couplings and 
using the intuition gained from Eq. (A.43). 

This model has the desired behavior, as it generates a PQ scale, fa = {S), of the order 
of the messenger scale. 


A.1.3 Type D Model 

From Eq. (A. 13) we have 


= 




'S21‘ 



(A.44) 


This model is equivalent to the one considered in Ref. [10] . The one-loop potential at order 

is 



(A.45) 
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where 


A = \S\^{Xl,+Xl,) + M^Xl^, (A.46) 

B = + 2A2^^m2|5P(A2^^ + A^^J + M^Xt,„ (A.47) 

C = A2^J5|^A2^^ + A^^J + M^\S\\Xl^ - 2XljXl^ + M^Xt^^, (A.48) 

D = |5|6(A2^^ - A2^J(A2^^ + A2^J2 - Xl^M^\S\^i3Xi, + 2Xl^Xl^ + SA^J 
+ 3Xt,,M^\S\\Xl^ - Xlj + Xl^M^. (A.49) 


At small I S'I 


vj.l>(S«M) = ^^ 


, , ,2 ^ 
1 + In A„,, 


+ 


327r2A2^^M4 


A 2 


A^,, - 4Af,, - Aj,, In 


^\31 

Xt,, M4 


+ 0 


M 6 


(A.50) 


The first non zero derivative is 


d^V^'^i\S\<€.M) 

95255*2 


Stt^XI^M^ 


- ^oxi. 


- K 22 111 


ALI^I 

A1,M4 


4\ n 


(A.51) 


This is positive for |5|/M —>■ 0, where the logarithmic term dominates. Thus we have a 
local minimum at the origin. At large |5| we have 


(|5| » M) = ... 2a2 ''mX2 ^ 
16 vr 2 (A 2 ^^+A 2 ^J 


F^Xt 


2A^22 + As2i 


I‘S'P(A^ 2 i + A^ ) 


A 2 


+A' 


S 22 






L 


+ AS 22 ) 


+ 0 


m 2 


(A.52) 


so the potential increases monotonically. We then infer that the minimum at the origin 
should be the global one. This can also be seen graphically in Fig. 4. 

We conclude, based on our analysis at one loop, that this model does not generate a 
large VEV for the field 5. The authors of Ref. [10] studied the same model at two loops 
and concluded instead that a large 5-VEV could be generated. 


B Anomalies and discrete symmetries 

In this Appendix we discuss the anomalies related to the DGS x Z^ introduced in 
Section 4 to address the issue of the axion quality. ‘ 

As we stated below Eq. (4.2) the orders M and N of the DGS must be chosen such 
that the expressions 

Z^^ : bqs = 0 (mod N ), (B.la) 

• aQx + bq§ = 2a + bqg = 0 {mod M), (B.lb) 

^Anomalies in the discrete symmetries used to protect the axion potential are considered in generality 
in non supersymmetric models in Ref. [30]. 
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Figure 4: A plot of the full Kahler effective potential (blue, Eq. (A. 45)) and the expansion 
of this potential around the origin (brown, Eq. (A. 50)) for model III. The choice of param¬ 
eters for this plot correspond to = -^521 = 1, Af = A = 1 and F = 

Note that here, for S/M > 0.2 the expansion ceases to be a good approximation to the full 
expression. 

are only satisfied for combinations of a and b such that a + 6 > 8. In the above expressions 
and label the charges of the chiral superfield under the and Z^ symmetries 
respectively. Without loss of generality we can choose the charge of the superspace coordi¬ 
nate q^ = 1; it follows that q^ = 2. Eor the sake of simplicity we take X uncharged under 
the Z^^ symmetry, qx = 0. Erom Eq. (B.l) a minimal choice of N and M, such that the 
constraint o -|- 6 > 8 is satisfied, would be A" = 3 and M = 5, with q^ = 1. 

Let us now consider the model dehned by Eq. (3.1). The form of the superpotential 
constrains the charges as follows 


qx 

+ 

q^i 

+ 

= 

= 0 (mod N ), 

(B.2a) 

qs 

+ 

q^i 

+ 

9^2 = 

= 0 (mod N ), 

(B.2b) 

qs 

+ 

q^2 

+ 


= 0 (mod N ), 

(B.2c) 

qx 

+ 

q^\ 

+ 

- 

y#i - 

= 2q^ (mod M), 

(B.2d) 

qs 

+ 

q^i 

+ 

^2 

= 2ql/ (mod M), 

(B.2e) 

qs 

+ 

q^2 

+ 

4, -- 

= 2q^ (mod M). 

(B.2f) 


Eor the anomalies we are interested in the charges of the fermions, which we denote with 
the subscript / next to the field subscript. Summing Eq. (B.2b) and Eq. (B.2c) we get 

i i 

while summing Eq. (B.2e) and Eq. (B.2f) we get 

-Qs- ‘^Qe ) = -^Qs (mod M). (B.4) 

i i 
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In these expressions the index i runs over all the messenger fields in the model. The 
requirement that our discrete symmetries be anomaly free amounts to ® 


where 


<l^iJ = 0 (“od r]N) , 

rf 

AG_a_^R = ^i{rf) q^.j + £ (adj G) = 0 (mod 7]m) , 

rf 


VN 


N for N odd , 
N/2 for N even , 


(B.5a) 

(B.5b) 


(B.6) 


The first sum runs over all irreducible representations t"/ of the gauge group G with Dynkin 
index £(rf), while the second sum runs over all fermions. In our conventions £{N) = 1/2 for 
the fundamental representation of SU(A^), and i (adj G) = N for the adjoint representation. 
For the PQ symmetry, using Eq. (B.3) we obtain 


^ (mod r/Tv) . (B.7) 

As we require qs 7 ^ 0, this is only satisfied for N even and qs = N/2. Thus we can 
choose a symmetry which is anomaly free. Then, to ensure a sufficient quality of 
the PQ symmetry without considering anomaly constraints on Z^j, the minimal choice is 
1^2^ X Zg, with qg = 2. 

For the /^-symmetry, using Eq. (B.4) we obtain 

^SU(5)-SU(5)-Z« = 5 — q§ = 0 (mod r]M) , (B.8) 

where we have included the contribution arising from the SU(5) gauginos. The only discrete 
i?-symmetry for M < 12 where this anomaly coefficient vanishes and all PQ violating non- 
renormalizable operators with a + 6 < 8 are forbidden is a symmetry where q^ = 5. 
These results, as stated in the main text, are dependent upon cancellations in only the 
R and PQ-breaking sectors resepectively. Namely, they will be different if one includes 
contributions from the visible sector or additional matter in the hidden sector. 
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